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A possible method to reconstruct the cosmic equation of state using strong gravitational 
\Q ■ lensing systems is proposed. The feasibility of the method is investigated by carrying out the 

reconstruction on the basis of a simple Monte-Carlo simulation. We show that the method 
can work and that the cosmic equation of state w(z) can be determined within errors of 
Aw ~ ±0.1 - ±0.2 when a sufficiently large number of lensing systems (N ~ 20) for z 1 
are precisely measured. Statistics of lensed sources in a wide and deep survey like the SDSS 
£C) ' are also briefly discussed. 

in ' 

■^j- ■ §1. Introduction 

o 

Because of recent progress in the capability of observational facilities, gravita- 
tional lensing phenomena in the high-redshift universe have come to play a very 
important role in the fields of cosmology and astrophysics. HP Gravitational lensing 
phenomena are useful not only to probe the nature of dark matter in the universe 
but also to test cosmological models. Recent reports on the measurement of the cos- 
mic shear field demonstrate this usefulness. l$'tl)'l§'El''B) It is also well recognized that 
strong lensing systems at high redshift are useful to constrain cosmological models, 
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in particular, to test the cosmological constant in the universe. I 

Recent observations of cosmic microwave background anisotropies and distant 
supernovae favor a spatially flat universe whose expansion is currently accelerating. 
The cosmological model that includes a cosmological constant is the simplest model 
that explains these observations. Motivated by those observational results, variants 
of the cosmological constant model have been proposed in the framework of the 
cold dark matter (CDM) cosmological model. For example, quintessence has been 
recently investigated. A model that includes quintessence has the attractive feature 
of explaining the 'coincidence problem' - thejieai. coincidence of the density of matter 
and the dark energy component at present. t!l)'[L3) 

Such generalized cosmological constants are regarded as a dark energy compo- 
nent in addition to the dark matter component in the universe. The dark energy 
component can be characterized by the equation of state w = Pq/pq, where pq is 
the pressure and pq is the energy density. If the dark energy is the cosmological 
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constant, then w = —1. For the quintessence model, w takes an arbitrary value sat- 
isfying — 1 < w < 0, and is a function of the redshift in general. Thus constraints on 
the equation of state are quite important for models of the dark energy component. 
Predictions of the CDM models with the dark energy co mponent n av e been tested 



by comparison with various cosmological observations.!^ >liJ)>[!j)>li3)>liZ|) An interest- 
ing approach to understanding the dark energy component is reconstruction of the 
cosmic equation of state or the quintessential potential. E3>B>E])>EiP 



Several authors havehwes 



jated gravitational lens systems as probes of the 
cosmic equation of state. c3"EjP'C3' Recently Futamase and Yoshida proposed a pos- 
sible method to probe the dark energy component using strong gravitational lensing 
systems. c3 They have pointed out that even one lensing system might give a strong 
constraint on the time variation of the dark energy component. Here, one might ask 
whether it is possible to reconstruct the cosmic equation of state or the quintessential 
potential using many strong lens systems, which may be detected using wide and 
deep surveys, e.g., the SDSS. In the present paper, motivated by the investigation 
in Ref. f24|), we investigate the feasibility of a scheme for reconstructing the redshift 
evolution of the cosmic equation of state. 

This paper is organized as follows. In §2, we briefly review theoretical formulas 
of the lensing model and a background cosmological model. In §3, a reconstruction 
scheme is worked out. In §4, the statistics of lensed sources detected in a wide 
and deep survey like the SDSS are briefly discussed. §5 is devoted to summary and 
conclusions. Throughout this paper we use units in which the velocity of light, c, 
equals 1. 

§2. Lens equation 

We restrict ourselves to a spatially flat FRW universe whose line element is 
written as 

ds 2 = a{ri) 2 {-dri 2 + dr 2 + r 2 dQ 2 {2) ), (1) 

where df2 2 2 ) is the line element of a unit two-sphere, rj is the conformal time and 
a(rj) is the scale factor normalized as a = 1 at present. We use a quintessential 
cosmological model consisting of a scalar field slowly rolling down its effect potential. 
Using the pressure pq and the energy density pq of the dark energy component, the 
effective equation of state of the dark energy component is written w = Pq/pq, 
which is a function of the redshift in general quintessential models. Assuming the 
equation of motion 

d( PQ a 3 ) + p Q d(a 3 ) = 0, (2) 
with the effective equation of state w(z), we have the solution 



PQ[Z) -(l + z) 3 exp 



} (z = 0) 



3 ' t&< 



m. (3) 



3 



In the case that the equation of state is written as w(z) = wo(l + z) u , we have 

'(l + zf-l 



f(z) = (l + z) d exp 



3w I 



(4) 



In the limit of a constant equation of state, v = 0, the dark energy density evolves 
as pq(z) oc f(z) = (1 + z J 3 0-+ w o)_ w e denote the density parameters of the matter 
component and the dark energy at present by J?o and &q(= 1 — A)), respectively. 
Then, from the Friedman equation, we have the formula for the angular diameter 
distance 



D A (zi,z 2 ) 



1 



dz' 



H (l + z 2 ) J Z1 [i7 (l + z') 3 + n Q f(z>)]V2 ' 



(5) 



where we assumed z\ < z 2 , and Hq =10Qh km/s/Mpc is the Hubble constant. 

Following previous investigations, IZP'ES we here consider the lens equation with 
a lens potential given by an isothermal ellipsoid model ascB 



$ = Ana: 



Dls 
D s 



1 - e)6\ + (1 + e)9l 



(6) 



where a v is the one-dimensional velocity dispersion, Dls is the angular diameter 
distance between lens and source objects, and Ds is the distance between the source 
and the observer. The ratio e of the minor axis to the major axis is related to the 
ellipticity e by e = + e)/(l — e). Thus the lens equation gives the elliptical 
image of the Einstein ring with minor axis and major axis as 



6± = 9 E VT±~e, 



(7) 



where 9e = ^Tra^DLS I Ds- Denoting the redshifts of the lens and the source by z\ 
and z s , respectively, we can write Dls = Da(zi, z s ) and Ds = Da(0, z s ). 



§3. Reconstruction scheme 



We assume a set of lensing systems whose observable quantities Z\ and z s and 
cr v , e and 0e are determined for each system. If a v , e and 0e are measured for 
a lensing system, Ds/Dls is determined from Eq. (|^). Thus the measurement of 
cr v , e and 6e determines the ratio of distances Ds/Dls- In our reconstruction 
procedure, we assume that a set of observable quantities zi, z s and R = Ds/Dls, 
instead of a v , e and 9e are determined. To be specific, we construct the set of N 
assumed lensing systems in a simulation as follows. First, we choose z s through a 
homogeneous random process in the range of redshifts zf 1111 < z s < z™ ax . Throughout 
this section we set z™ m = 0.5 and z™ ax = 1. In a similar way, we choose z\ in the 
range z™ in <z\< zf 1 ^, where we set zf 1111 = 0.1 and z ; max = z™ ax — 0.2 in the present 
paper. Using a cosmological model, we then compute Ds/Dls using the distance- 
formula (5). Observational errors will cause an observed value of Ds/Dls to differ 
from the theoretical prediction. We assume that the ellipticity and the angle of the 
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Einstein ring are well determined for each lensing system and that the observational 
noise is dominated by an error in a v . Then we set 



R 



Ds(zi,z s ) 
Dls(zi,z s 



x (1+p) 2 , 



(8) 



where p is a random variable that is characterized by a Gaussian probability function 
with mean and the variance a 2 . □ 

We next consider reconstructing the cosmic equation of state from the set of N 
lensing systems specified by zn, z s i and Ri, where the subscript i denotes the data 
of the i-th lensing system. Our method is a simple application of the reconstruction 
scheme from supernovae data proposed by Chiba and Nakamura.cil We adopt the 
following fitting formula for the angular diameter distance: 



Da{(zi,z 2 ) 



1 



1 



H l + z 2 



m(zi) - vdz2. 



with 



m{z> 



2a[y- 8 + (3y- 6 + 1V ~ A + 5y~ 2 + C]- 1/8 , 



(9) 



(10) 



and y = 1/Vl + z. Here the subscript f refers to fitting formulas. We use the 
constraint equation from the Friedman equation© 



a 



[1 + f3 + 7 + S + C] 9/8 



1 + 3/3/4 + 7 /2 + 5/4 
Then we determine the parameters j3, 7, 5 and £ so as to minimize 

(Ri - Rfi(zn, z si , P, 7, ^ C)) 2 



(11) 



N 

E 

i=l 



F A (P,j,6,0 
where Ru is defined by 

Rtt(zu, z si , {3,7,6, () 



R 2 



(12) 



DAf(0,z s 



??f(o) - m(zsi) 



D Ai (zu, z si ) Vf(zn) ~ m(z S i) ' 



(13) 



Thus Ra does not depend on a, so the relation (|ll]) is needed to determine a in our 
prescription. Finally, the cosmic equation of state is given by 



Wf(z) 



Ayd 2 r] { /dy 2 



3(dr H /dy)[f2 (dr 1i /dyy 



(14) 



Figure 1 demonstrates the feasibility of the reconstruction scheme for cosmolog- 
ical models whose parameters are Qq = 0.3, wq = —0.5, v = 0.4 (left panels) and 
i?o = 0.3, wq = —0.7, v = —0.4 (right panels). In each panel the dashed curve 



* Note that the variance a 2 is distinct from a 2 ,, which represents the velocity dispersion for the 
lens model. 
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represents the theoretical curve w(z) = wq(1 + z) v , while the solid curve and shaded 
region represent the mean and the 1-sigma variance of the ensemble average of the 
reconstructed curves. The upper panels display the results for which the parameters 
(3, 7,5 and £ are determined by minimizing the function Fa- The middle panels 
correspond to the case in which the parameters are determined so as to minimize 
the function Fb instead of Fa , where 



The lower panels display the results for which the reconstructed curve of Wf(z) is 
obtained by averaging the two reconstructed curves obtained using Fa and Fb- In 
Fig. 1 we assumed a = 0.03 and N = 20. 

The capability of the reconstruction scheme depends on the errors involved in 
measuring the lensing systems, which are incorporated by p (or a) in Eq. (g) in our 
simulation. Figure 2 displays the variance of reconstructed equation of state Aw{z) 
at z = 0.5 as a function of a for the same cosmological models as in Fig. 1, i.e., with 
i?o = 0.3, wq = —0.5, v = 0.4 (left panels) and i?o = 0.3, wq = —0.7, v = —0.4 
(right panels). This shows that the equation of state can be reconstructed within 
errors of Aw ~ ±0.1 - ±0.2. 

In our simulation, for a set of data obtained through a random process, the 
scheme fails in reconstruction and gives a non-realistic value of Wf(z). We removed 
such a case in our simulation and considered only the case that the reconstructed 
equation of state satisfies —3/2 < Wf(z) < for < z < 1. 'Failed' reconstructions 
appear more frequently as a becomes larger than 0.03. For example the 'failed' 
reconstruction appears at the rate of 3%, 20% and 50% for a = 0.01, 0.03 and 0.05, 
respectively. Hence our reconstruction scheme might not be feasible for a £ 0.03 for 
general models other than the quintessential model. 



In the previous section we have assumed a set of N (= 20) lensing systems 
distributed in the range 0.5 < z s < 1. In this section we give some gromids for 
this assumption. Lensing statistics have been thoroughly investigated,!^ EZP and 
following previous investigations, we can roughly estimate the number of lensed 
sources per redshift as (see the Appendix for a derivation) 



v 



Rfi(zu,Zsi,0,l,8, 2 



F B (J3,i,6,C) = Y, 



(15) 



§4. Discussion — Lensing statistics 



dN(z a ) 
dz s 



5.8 x 10 7 (H r(z s )) 



5 dH r(z s ) n ga \(z s ) 
dz s h 3 Mpc" 



-3 




(16) 



where r(z) is the comoving distance at redshift z, n ga \(z s ) denotes the comoving 
number density of galaxies, and nhalo denotes the comoving number density of the 
lensing halo, which is assumed to be constant throughout the universe. In deriving 
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Eq. (16), the singular isothermal sphere model is assumed as the lens model, and 
the probability that multiple images appear is computed. 

Figure 3 displays the expected distribution of lensed sources dN/dz s for cosmo- 
logical models with Qq = 0.3, wq = —0.5, v = 0.4 (left panels) and Qq = 0.3, wo = 
—0.7, v = —0.4 (right panels). Here we adopted a luminosity function of the APM 
galaxies E9) fitted to the Schechter function 

4>{L)dL = <p(jL\ expf-^U^, (17) 

with (f>* = 1.40 x 10~ 2 /i 3 Mpc -3 , k = -0.97, and M* = -19.50 + 51og 10 /i. Then 
the comoving number density of galaxies at z which are brighter than the limiting 
magnitude B\ im is given by 

n gal (z,<B hm ) = / <j>{L)dL = <t>*r[K+l,x{B lim ,z)], (18) 

JL(B Um ,z) 

where x(B hm , z) = L{B hm ,z)/L* = [d^/lh- 1 Mpc] 2 10 2 - 2 -°- 4Bli ™, r[n,x] is the 
incomplete Gamma function, and di{z) is the luminosity distance. In each panel of 
Fig. 3, the three solid curves correspond to the cases By im = 21, 22, and 22.5. I 

The expected number of lensed objects N is obtained by integrating dN/dz s . 
For example we have N ~ 410, 2800, 7000 for B hm = 21, 22, 22.5, respectively, 
for the model with Qq = 0.3, wq = —0.7, v = —0.4 within a solid angle of ir str. 
This is almost the same result as that for the model with = 0.3, wq = —0.5 and 
v = 0.4. We thus conclude the SDSS project in progress should detect many lensed 
objects. Of course, the above is an over-estimation because the formula contains the 
case in which the ratio of the magnitude of the two lensed images is infinite. The 
actual number of 'clean' lensed systems near the Einstein ring should be significantly 
smaller. If 1% of the above estimation are clean lensed systems, this number should 
be sufficient to carry out the reconstruction proposed in this paper. 

§5. Conclusion 

In this paper we have proposed a possible method to reconstruct the cosmic 
equation of state using strong gravitational lensing systems. We have investigated 
the feasibility by working out a reconstruction process based on a simple Monte- 
Carlo simulation. Our result shows that the method can work when a sufficiently 
large numbers of lensing systems (N ~ 20 at z ^ 1) are precisely observed. For 
example, if the velocity dispersion (the lens model) can be measured within errors 
smaller than ~ 3%, the cosmic equation of state can be determined within errors 
Aw ~ ±0.1 - ±0.2. The accuracy required for the measurement of the velocity 
dispersion might be quite high. However, it should be noted thaLsuch accurate 
measurements have been made using by the Keck-II 10m telescope. IH£P'E2P 

* In the present paper we neglect the magnification bias, because galaxies in the range of rather 
small redshift, 0.5 < z 3 < 1, are considered as lensed sources. The magnification bias enhances 
the expected number of the lensed sources. However, this effect does not alter our conclusion 
qualitatively. 
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We have adopted a simple model for the redshift distribution function of lensed 
sources. Specifically, in the present paper we have assumed a homogeneous random 
distribution in the range 0.5 < z s < 1. The feasibility of the reconstruction scheme 
sensitively depends on the redshift distribution of the lensing system. Hence, further 
investigations of distributions of lensing systems are needed. In general, the capa- 
bility of our reconstruction scheme becomes worse as the redshift of lensed sources 
becomes higher. This feature originates from the fact that the fraction of the dark 
energy component becomes smaller relative to the matter component as the redshift 
becomes higher. This makes Da(zi, z s ) less sensitive to the cosmic equation of state 
w(z) as zi and z s become larger. 

Our investigation is based on some other idealistic assumptions. We have as- 
sumed a flat spatial geometry of the universe. The crucial assumption might be that 
the matter fraction in the universe Q$ is determined precisely. The errors in the 
determination of i?o cause additional errors in reconstructing w(z). 

We have assumed a singular isothermal ellipsoid as a model of our lensing galax- 
ies because of its simplicity. However, it is not known if the singular isothermal 
ellipsoid model is a good model for lensing galaxies. In fact, the modeling of lensing 
galaxies is not easy, because of many theoretical and observational ambiguities such 
as the existence of an anisotropic velocity dispersion. However, there has been some 
progress in this direction. Dynamical observations and models of lensing elliptical 
galaxies in the local universe have been studied by Rix et al.H' Furthermore, the 
relation between the observed stellar velocity dispersion and the velocity dispersion 
associated with dark matter has been studied in the case of singular isothermal 
sphere model.© It is expected that such studies will provide more accurate results 
in the near future, based on the large sample of galaxies in SDSS. Thus we hope that 
our method will arrow for realistic reconstruction of the cosmic equation of state in 
the near future. 
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Here we briefly review the estimation of lensing statistics. E^P'lHZP Introducing the 
cross section na^ T , the optical depth due to the gravitational lensing can be written 
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Appendix A 
Lensing Statistics 




(A-1) 



s 



where nhaio(-zz) is the comoving number density of the halo at redshift z\. Adopting 
a singular isothermal sphere as the model of the lens halo, we have 



where r(z) is the comoving distance and a v is the velocity dispersion. In deriving 
Eq. ( |A-2j ) , we have used the assumption of a spatially flat universe and the fact that 
a cr is written as 

a CI = (1 + z i )- 1 r(^)47ra2^, (A-3) 

J->s 

where a cr is the critical radius for multiple images. 

The probability function that a source at z s becomes multiple images, which is 
defined by P (z s ) = Jq" dr, can be written as 

Pfe, = 2.2 X 10-3( gor (,.)) a ( 10 _ 2 ^ pc _ 3 ) (^) 4 . (A-4) 

where we have assumed the comoving number density is constant n^^z) = nhaio- 
Then the expected number of lensed galaxies per unit solid angle is computed by 

N= r dz s n g&l (z s )P(z s )r(z s ) 2 ^p±, (A-5) 
Jo dz s 

which reduces to 

,7 T, (l j , x , 5 dH r(z s ) n ga iO s ) 



2 



poo 

N = 5.8 x 10 7 / dz s (H r(z s )f ., _. 

jo dz s h 6 Mpc ' 

\ / \ 4 

(A-6) 



10-2/i 3 Mpc" 3 / V 250 km/s 
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Fig. 1. Feasibility of the reconstruction scheme for cosmological models whose parameters are J?o = 
0.3, wo = —0.5, v = 0.4 (left panels) and Oo = 0.3, wo = —0.7, v = —0.4 (right panels). Here 
we have assumed lensed sources (N = 20) distributed in the range 0.5 < z s < 1 and a = 0.03. 
In each panel, the dashed curve is the theoretical curve, while the solid curve and shaded region 
represent the mean and the 1-sigma variance of the reconstructed curves. The upper (middle) 
panels display the results for which the parameters /3, 7, 8 and £ are determined by minimizing 
the function Fa (Fb)- The lower panels display the result for which the reconstructed curve 
Wf(z) is obtained by averaging the reconstructed curves obtained from both Fa and Fb- 
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. 2. Variance of the reconstructed cosmic equation of state Aw at z = 0.5 for cosmological 
models whose parameters are f2o = 0.3, wq = —0.5, v = 0.4 (left panels) and Qo = 0.3, wq = 
—0.7, v — —0.4 (right panels). Here we have assumed that the number of lensing sources N = 20 
is distributed in the range 0.5 < z s < 1, as in Fig. 1. 
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Fig. 3. dN/dz s for cosmological models with Qq = 0.3, too = —0.5, v = 0.4 (left panels) and 
(2o = 0.3, too = —0.7, v = —0.4 (right panels). Here we have adopted the luminosity function 
for galaxies fitted to the Schechter function. In each panel, the solid curves correspond to the 
cases in which the limiting magnitude is Bn m = 21, 22 and 22.5, from bottom to top. 



